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The  Fractal  Nature,  Graph  Invariants,  and  Physicochemical 
Properties  of  Normal  Alkanes 
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department  of  Physics  and  Atmospheric  Sciences,  Jackson  State  University, 
Jackson,  Mississippi  39217 

Abstract 

The  variation  in  the  physicochemical  properties  of  linear  normal  alkane 
molecules  with  the  graph  invariants  carbon  number  and  the  Wiener  index  is 
investigated.  It  is  demonstrated  that  the  variation  can  be  interpreted  in  terms 
of  the  fractal  dimension  of  these  species.  A  general  formula  expressing  the 
precise  nature  of  the  relationship  has  been  derived. 

Introduction 

As  the  most  extensively  studied  of  all  the  homologous  series,  the  alkanes 
may  now  be  said  to  be  generally  well  understood.  Extensive  compilations  of 
data  on  their  physicochemical  properties  have  been  amassed^  and  numerous 
structure-property  correlations  have  been  performed  using  a  wide  variety  of 
different  graph  invariants. 2  Such  correlations  tend  to  have  only  a  limited  range 
of  applicability,  however,  and  attempts  to  extend  the  range  usually  lead  to 
excessively  complicated  expressions. 3  One  important  reason  for  the  difficulties 
encountered  in  modelling  the  alkanes  in  this  way  is  that,  for  any  given 
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physicochemical  property,  the  increments  to  that  property  from  each  subsequent 
addition  of  a  methylene  (CH2)  unit  along  the  series  are  not  constant.  The 
contributions  from  the  methylene  components  become  progressively  smaller 
as  the  carbon  number,  n,  increases.^  This  type  of  behavior  is  evident  from  Figure 
1,  which  shows  the  plot  of  boiling  point  against  £  for  normal  alkanes  up  to  the 
C40  member. 

In  recent  years  many  different  graph  invariants  have  been  proposed  for  the 
characterization  of  chemical  species.  These  are  usually  referred  to  in  the  chemical 
literature  as  topological  indices.®'®  Such  indices  reflect  in  different  ways  the 
size  and  shape  of  the  molecules  they  characterize  and  also  provide  some  measure 
of  the  degree  of  branching  present.^  The  relationship  between  such  indices  and 
the  physicochemical  properties  of  the  alkanes  is  clearly  not  a  linear  one  (note 
the  pronounced  curvature  in  Fig.  1).  Even  if  the  plots  are  made  on  logarithmic 
scales,  some  of  the  curvature  remains,  as  will  be  apparent  from  Figure  2,  which 
shows  the  variation  of  In  £  with  In  (boiling  point)  for  normal  alkanes.  To  explain 
this  curvature  we  shall  make  use  of  fractal  geometry  and  demonstrate  that  this 
curvature  arises  from  the  fractal  nature  of  the  alkanes  themselves. 

The  concept  of  the  random  fractal  has  been  used  to  characterize  any  object 
whose  form  is  extremely  irregular,  ramified  and/or  fragmented  at  length  scales 
on  which  a  self-similarity  exists.®  Widespread  use  of  the  concept  has  been  made 
in  both  chemistry  and  physics.^" ^  Although  the  full  description  of  a  fractal 
involves  several  different  universal  scaling  functions  and  exponents,  such  as 
the  fractal  dimensionality,  the  topological  dimensionality,  the  order  of 
ramification,  the  connectivity,  and  lacunarity”'®,  we  shall  focus  here  only  on 
the  fractal  dimensionality.  The  fractal  dimensionality  is  one  of  the  most  important 
scaling  exponents  and  quantifies  the  manner  in  which  the  mass  of  an  object 
increases  with  its  size.**^  If  m  is  the  mass  of  the  object  and  £  its  radius  of  gyration. 


we  may  write: 


The  carbon  number,  n,  in  an  alkane  molecule  was  first  employed  over  100 
years  ago  as  a  means  of  characterizing  these  molecules.15  it  represents  the 
simplest  possible  graph  invariant,  for  it  may  be  defined  as  the  number  of  vertices 
in  the  hydrogen-suppressed  graphs  of  alkane  species.^  As  such,  the  carbon  number 
may  be  deemed  to  be  a  topological  index,  though  it  is  one  which  has  very  low 
discriminating  power.  All  isomers  corresponding  to  a  given  value  of  ii  will  have 
the  same  index.  Such  an  index  is  thus  completely  unsuitable  for  the 
characterization  of  branched  species.  Since  this  index  came  into  vogue,  many 
other  topological  indices  have  been  developed  with  the  specific  aim  of 
characterizing  the  branching  in  molecular  species.^  We  now  describe  one  such 
index  which  we  shall  use  in  our  investigation  of  the  behavior  of  linear  alkanes. 

The  topological  index  we  employ  was  introduced  by  Wiener*^  and  is  nowadays 
commonly  referred  to  as  the  Wiener  index.  Although  the  index  was  first  put 
forward  in  1947,  its  manifold  applications  in  chemistry,  physics,  and  several 
other  sciences  have  only  recently  come  to  the  fore.  Details  of  the  index  and 
its  applications  are  to  be  found  in  a  review  by  one  the  present  authors.*^  The 
Wiener  index  may  be  defined  as  one  half  the  sum  of  the  entries  in  the  distance 
matrix,  D(G),  of  the  graph  G  of  the  molecule  under  consideration.  Normally 
C  is  taken  as  the  hydrogen-suppressed  graph  of  the  species.  Thus,  as  will  be 


apparent  from  Figure  3,  the  Wiener  index,  W(G),  for  the  molecule  of  normal 
butane  is  10.  For  normal,  i.e.  unbranched,  alkanes,  Wiener  proved1*’  the  following 
general,  closed  formula  for  W(G): 


W(G)  =  -g  (n3  -  n) .  (2) 

We  shall  concern  ourselves  here  only  with  normal  alkane  species.  In  Figures 
4  and  5  plots  of  W(G)  versus  the  boiling  points  of  normal  alkanes  up  to  C40  member 
are  presented  using  respectively  ordinary  and  logarithmic  scales. 

Modelling  the  Curvature 

The  pronounced  curvature  evident  in  Figures  1  and  4  is  typical  of  plots  of 
most  physicochemical  parameters  against  some  appropriate  graph  invariant  used 
to  characterize  the  structure  of  species.  An  attempt  to  model  relationships 
of  this  general  type,  led  Walker1**  to  postulate  a  biparametric  equation  of  the 
following  form: 

P  -  a[W(C)lb  (3) 

where  P  is  some  physicochemical  property,  and  a,  and  J3  are  constants  to  be 
determined  for  the  range  of  values  of  the  property  P  under  consideration.  For 
limited  ranges,  estimates  of  a  and  kj  are  commonly  obtained  from  logarithmic 
plots  such  as  that  shown  in  Figure  5.  For  more  extended  ranges,  however,  it 
is  clear  that  this  method  is  not  applicable,  for  the  plots  are  then  far  from  being 
straight  lines. 

Since  the  first  few  points  of  the  plot  in  Figure  5  (corresponding  roughly  to 
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1  <  £  <  10)  may  be  fitted  reasonably  well  on  one  straight  line,  and  the  last  few 
points  (corresponding  to  30  <  n  <  40)  may  be  fitted  on  another  straight  line,  it 
is  of  interest  to  inquire  if  the  two  extreme  slopes  have  physical  significance. 
Let  the  slope  of  the  first  few  points  be  ^-|  and  that  of  the  second  few  points 
be  b2.  Clearly  we  have  b-]  >  and  the  range  of  b  values  lying  between  b-j 
and  b2  will  correspond  to  a  crossover  region  characterized  by  a  transition  from 
one  power  law  behavior  (having  exponent  b-j)  to  another  power  law  behavior  (having 
exponent  bj).  The  Walker-type  relationship  in  equation  (3)  should  thus  more 
correctly  be  written  as  several  different  relationships.  For  the  extreme  cases 
referred  to  above,  the  two  relationships  may  be  written  in  the  form: 

P  =  a-|  [W(G)^1  (n  small)  (4a) 

P  =*  £2  [W(G)£2  (£  large)  (4b) 

where  a-j  and  £2  are  appropriate  constants. 

The  type  of  behavior  described  by  equations  (4a-4b)  has  previously  been 
observed  and  commented  upon  by  several  authors.^-21  The  exponents  Jd-j  and 
b?,  which  we  shall  focus  on  exclusively  hereafter,  can  be  interpreted  in  terms 
of  the  fractal  nature  of  normal  alkane  species.  Detailed  studies  of  the  alkanes 
near  to  their  boiling  point  have  led  to  the  conclusion  that  the  molecules  behave 
more  or  less  as  rigid  rods  when  they  are  short22  and  as  random  chains  when  they 
are  very  long. 23  For  our  purposes  it  is  thus  reasonable  to  assume,  as  a  good 
first  approximation,  that  short  normal  alkane  molecules  (1  <  £<  10)  are  completely 
stiff,  whereas  long  normal  alkane  molecules  (jn  >  35)  are  completely  flexible. 
It  is,  of  course,  well-known24  that  £  =  35  will  not  be  long  enough  to  yield  a 
completely  flexible  chain;  this  figure  was  selected  somewhat  arbitrarily  because 


of  the  paucity  of  boiling  point  data^  on  species  having  £  >  40.  In  fact,  ri  =  30 
lies  in  the  crossover  region  between  rigid  rod  and  flexible  chain  type  behavior. 24 
Strictly  speaking,  therefore,  precise  behavior  of  one  type  or  the  other  cannot 
be  ascribed  to  alkanes  in  this  region. 

Alkane  Fractal  Dimensionalities 

In  our  model  we  thus  imagine  short  alkane  molecules  to  exist  in  the  form 
of  stiff  rods,  as  illustrated  in  Figure  6.  The  end-to-end  length,  £p,  of  such  a 
chain  has  been  shown^S  to  be  proportional  to  n: 

In  °  H  •  (5) 

On  the  other  hand,  the  flexible  chain  characteristic  of  long  alkanes,  which  we 
illustrate  in  Figure  7,  is  known^4  to  have  an  end-to-end  length,  £p,  that  assumes 
the  following  proportionality: 

rp  a  ni  .  (6) 

There  is  thus  a  clearcut  change  in  the  power  law  which  obtains  in  going  from 
a  stiff,  rodlike  short  alkane  chain  to  a  flexible,  gaussian,  random-walk-like  long 
alkane  chain. 

From  equation  (1),  we  may  therefore  conclude  that  for  small  values  of  ji, 
df  =  1,  whereas  for  large  values  of  £,  df  =  2.  The  difference  in  these  fractal 
dimensionalities  reveals  that  we  are  dealing  with  two  self-similar  configurations 
of  the  alkanes  based  on  two  different  topological  length  scales.  These  are  (i) 
a  rodlike  configuration  for  small  £  which  persists  to  some  upper  cutoff  value 
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ry;  and  (ii)  a  flexible  chainlike  configuration  for  large  _n  where  n  >  _nc.  It  is  these 
two  different  fractal  dimensionalities  which  we  use  to  explain  the  curvature 
observed  in  the  plots  of  physicochemical  properties  of  the  alkanes  versus  the 
carbon  number,  n>  or  indices,  such  as  the  Wiener  index,  which  correlate  closely2^ 
with  n.  As  is  evident  from  equation  (2),  the  Wiener  index  for  large  £  may  be 
expressed  to  a  good  approximation  as  follows: 

W(G)  a  <n3  -  n)  =  n3  (1  -  rf2)  -  n3  .  (7) 

Rearrangement  of  the  two  expressions  given  in  equations  (4a-4b)  to  model 
the  variation  of  the  physicochemical  property,  P,  with  the  Wiener  index,  yields 
the  new  equations: 

2 

W(G)  a  P-1  (n  small)  (8a) 

and 

2 

W(G)  a  P^2  (n  large).  (8b) 

Because  the  physicochemical  properties  of  the  alkanes  will  depend  in  general 
upon  the  configuration  of  the  molecules  and  the  configuration  adopted  will  in 
turn  depend  upon  n,  P  can  be  expressed  in  terms  of  the  following  relationship: 

P  a  nil  ,  (9) 

where  k_  is  an  exponent  which  we  suppose  assumes  the  values  k_1  and  k_2  thus: 
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JiVjil 

W(G)  an  (n  small) 


(10a) 


k  ?/b2 

W(G)  an  (n  large) 


(10b) 


If  no  cognizance  is  taken  of  the  actual  configurations  adopted  by  the  alkane 
species,  i.e.  if  no  account  is  taken  of  the  differing  self-similarities  and  their 
fractal  dimensionalities,  it  is  possible  to  derive  two  further  relationships.  Using 
the  value  of  W(G)  given  in  equation  (7),  and  comparing  it  with  the  exponents 
in  equations  (10a- 10b),  we  obtain  the  relationships: 


kj/Jbi  =  3 


(n  small) 


(11a) 


k2/b2  =  3 


(£  large)  . 


(11b) 


To  relate  the  results  in  equations  (11a-11b)  to  the  fractal  dimensionality,  df, 
of  rp,  we  make  use  at  this  point  of  an  appropriately  modified  form  of  the 
expression  in  equation  (1): 


W(G) 


bi/k_i 


W(G) 


b?/k? 


(n  small) 


(n  large) 


(12a) 


(12b) 


Scaling  Theory  Applied  to  the  Alkanes 


To  interrelate  the  two  sets  of  exponents  {b-j,  k^}  and  {  bj,  k£},  which  we 
have  defined  as  being  valid  in  the  two  limiting  cases  (i)  n  <  nc  and  (ii)  £  >  £(-, 
use  is  now  made  of  scaling  arguments  similar  to  those  presented  by  de  Gennes.27 
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The  size  of  the  molecules,  characterized  by  r n,  will  clearly  be  different  for  the 
two  cases.  Only  at  the  point  where  =  ric  will  the  two  cases  by  precisely  equal, 
namely  when 

In  ( n. <  £ic)  =  £p  (n_  >  He  •  (13) 


At  this  point,  it  then  follows  from  equations  (12a-12b)  that: 


kVJil^fl  b?/k9df? 

W(G)  “  =  W(G) 


and  therefore  that: 


kl  =  J 22 

ill  dfi  k 2  df 2 


(IS) 


Since  from  our  earlier  discussion  we  already  know  that  df-j  *  1  and  df2  =  2,  we 
are  now  in  a  position  to  state  our  final  new  result  as: 


kl 

ii 

b2 

2k2 

irimental  Test  for  the  Derived  Relationship 


(16) 


To  test  whether  the  new  relationship  derived  above  in  equation  (16)  is  a  valid 
one,  determinations  were  made  of  the  various  slopes  exhibited  in  the  plot  shown 
in  Figure  5.  Using  a  least  squares  regression  analysis,  we  computed  the  slope 
of  the  first  nine  points  in  the  Figure  (corresponding  to  C2  through  Ciq  alkanes), 
the  slope  of  the  last  nine  points  (corresponding  to  alkanes  C32  through  C40), 
as  well  as  the  slopes  of  several  other  sets  of  nine  points  chosen  from  the  Figure. 


-10- 


The  slope  for  the  first  set  of  nine  points  was  found  to  be  0. 1133°K,  and  that 
for  the  last  set  equalled  0.1753°K.  If  each  of  the  slopes  obtained  is  divided  by 
the  slope  for  the  first  set  of  points  (0.1133°K),  the  ratios  should  converge  to 
the  value  0.S  as  £  increases.  A  plot  of  these  ratios  against  the  average  value 
of  ri  represented  by  each  set  of  nine  points  chosen  is  presented  in  Figure  8. 

Unfortunately,  reliable  boiling  point  data  for  the  alkanes  are  not  available 
beyond  the  C4Q  member,  and  so  this  perforce  represents  an  abrupt  cutoff  point 
for  our  plot  in  Figure  8.  From  this  plot,  however,  it  may  be  seen  that  whereas 
the  ratio  byb2  starts  out  from  the  precise  value  of  unity,  by  the  time  the  last 
set  of  nine  points  is  reached  (corresponding  to  an  average  C35)  the  ratio  has 
fallen  to  the  value  0.646.  It  seems  highly  likely  that  this  fall  toward  the 
theoretical  limit  of  0.5  would  continue  if  sufficient  boiling  point  data  were 
available  to  extend  our  plot  beyond  this  artificially  imposed  limit. 

Conclusion 

A  formula  has  been  derived  showing  how  the  fractal  dimensionalities,  df, 
for  £p  characterizing  normal  alkane  species  are  related  to  the  exponents  b-]  and 
b?  in  Walker-type  plots  of  In  P  versus  In  W(G)  for  ri  small  and  £  large.  The  entirely 
reasonable  assumption  has  been  made  that  the  two  extreme  df  values  are  one 
for  short-chain  normal  alkanes  and  two  for  long-chain  normal  alkanes.  We  have 
demonstrated  that  the  ratios  of  the  slopes  obtained  from  plots  in  In  P  versus 
In  W(G)  should  asymptotically  approach  the  theoretical  limit  value  of  0.5  as  jn 
becomes  very  large.  Analysis  of  these  plots  for  normal  alkanes  up  to  and  including 
the  C40  member  reveal  that  the  ratio  falls  in  value  from  unity  for  the  C2-C10 
sets  of  points  to  0.646  for  the  C32-C40  set  of  points.  It  appears  very  likely  that 
the  observed  fall  in  value  of  the  ratio  would  continue  if  data  were  available  for 
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chains  containing  more  than  forty  carbon  atoms,  so  that  the  theoretical  value 
of  0.5  would  be  more  closely  approached.  From  our  results  we  may  therefore 
conclude  that  the  Wiener  topological  index  is  able  to  provide  valuable  information 
on  the  state  of  randomness  of  normal  alkane  chains  at  the  boiling  point.  The 
general  approach  we  have  outlined  here  could  clearly  be  extended  to  the  study 
of  a  variety  of  other  systems  and  many  different  physicochemical  parameters. 
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Figure  1.  Plot  of  boiling  point  temperature  against  carbon  number  for  normal 
alkanes  from  the  C-]  to  the  C40  member. 

Figure  2.  Plot  of  natural  logarithm  of  boiling  point  temperature  against  natural 
logarithm  of  the  carbon  number  for  normal  alkanes  from  the  C-]  to 
the  C40  member. 

Figure  3.  Schematic  derivation  of  the  Wiener  index  for  the  molecule  of  normal 
butane. 


Figure  4.  Plot  of  the  boiling  point  temperature  against  Wiener  index  for  normal 
alkanes  from  the  C2  to  the  C40  member. 

Figure  5.  Plot  of  the  natural  logarithm  of  boiling  point  temperature  against 
natural  logarithm  of  the  Wiener  index  for  normal  alkane  species  from 
the  C2  to  the  C40  member. 

Figure  6.  An  illustration  of  the  end-to-end  length,  rp,  for  short,  stiff,  normal 
alkane  species  (1  <  £  <  10). 

Figure  7.  An  illustration  of  the  end-to-end  length,  rp,  for  long,  flexible,  normal 
alkane  species  (£  >  35). 

Figure  8.  Plot  of  the  slopes  for  different  sets  of  nine  points  (taken  from  Figure 
5)  against  the  average  value  of  £  represented  by  each  set  of  points. 
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